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Guiding-center polarization and magnetization effects in gyrokinetic theory 

Alain J. Brizard 

Department of Physics, Saint Michael's College, Colchester, VT 05439, USA 

Higher-order guiding-center polarization and magnetization effects are introduced in gyrokinetic 
theory by keeping first-order terms in background magnetic-field nonuniformity. These results con- 
firm the consistency of the two-step perturbation analysis used in modern gyrokinetic theory. 

o ' 

■ I. INTRODUCTION 

\ Polarization and magnetization effects play a fundamental role in modern gyrokinetic theory [J, Q ■ The standard 
form of modern gyrokinetic theory is derived by a two-step Lie-transform perturbation analysis that retains the 
effects of first-order (eg) guiding-center drifts associated with the nonuniformity of the background magnetic field and 
gyrocenter effects up to second-order (e^) in electromagnetic- field fluctuations that perturb the background guiding- 
i— |" center plasma. First-order (e^) gyrocenter polarization and magnetization effects (which result from terms of order ej 
tIh , in the gyrocenter Hamiltonian) are fully retained in nonlinear gyrokinetic theory. Because the ordering parameters cb 
Oh- and eg are often comparable (eg ~ eg) in many practical applications of gyrokinetic theory, however, it is sometimes 
argued 3] that first-order (es) guiding-center polarization and magnetization effects (which result from terms of order 
€b £(5 in the gyrocenter Hamiltonian) should also be retained in nonlinear gyrokinetic theory for a consistent treatment 
5j ' of polarization and magnetization effects. 

The two-step derivation of modern gyrokinetic theory is based on a sequence of two near-identity phase-space 
transformations z — > z — > z from local particle coordinates z = (x, Vq\\, fJ,Q,6o) to guiding-center coordinates z = 
(XjW|,/x, 0) and then to gyrocenter coordinates z = (X, U||,p, &). The purpose of the guiding-center transformation 
4, 5] zo — >• z (with small parameter eg) is to asymptotically decouple the fast gyromotion of charged particles in a 
strong weakly- nonuniform (eg <C 1) background magnetic field and construct the guiding-center magnetic moment 
fi = /j,q + es fii + ■ ■ ■ as an adiabatic invariant (from the local particle magnetic moment fig = m\vo±\ 2 /2B). The 
I introduction of low-frequency electromagnetic-field fluctuations destroy the adiabatic invariance of the guiding-center 
l— "™ '■ magnetic moment fj,, which requires the gyrocenter transformation z — > z (with small parameter e^) in order to restore 
I \ the adiabatic invariance of the gyrocenter magnetic moment ~p = /i + eg fJ 1 + ■ ■ ■ . 

■ Each transformation introduces a polarization charge density in the gyrokinetic Poisson equation and polarization 
and magnetization current densities in the gyrokinetic Ampere equation Q. These effects explicitly involve the 
generalized gyroradius p = x — X defined as the displacement of the gyrocenter position X from the particle position 

OO ■ x. As a result of the guiding-center and gyrocenter transformations, the generalized gyroradius 

cn ; 

P = (Pogc + cb Pi gc + • • • ) + (eg Pi gy + • • • ) (1) 

is decomposed into the guiding-center gyroradius p gc = p 0gc + es Pi gc + • • • and the gyrocenter gyroradius 
p gy = eg Pigy + ••• (which vanishes in the absence of fluctuations). As will be discussed below, polarization ef- 
fects are associated with the gyrorangle-averaged displacement (p) while magnetization effects arc associated with 
k>( , the gyrorangle-dependent displacement p = p — (p) . 

The first-order (eg) gyrocenter polarization and magnetization effects associated with p lgy (which can also include 
terms of arbitrary order in e^) have been discussed elsewhere (HQ. Moreover, we note that, since electromagnetic- 
field fluctuations satisfy the gyrokinetic ordering |p 0gc • V| ~ 1 (i.e., perpendicular wavelengths are of the same order 
as the lowest-order gyroradius), the guiding-center gyroradius-expansions considered in the present work apply only 
to the derivations of guiding-center polarization and magnetization effects. Hence, the first-order (eg) gyrocenter 
polarization and magnetization effects will be ignored in what follows. 

Because the lowest-order guiding-center gyroradius is explicitly gyroangle-dependent (i.e., (Po gc ) = 0)i there is 
no zeroth-order term in the guiding-center polarization while the guiding-center magnetization has a non-vanishing 
zeroth-order term. First-order guiding-center polarization effects were investigated previously |7H10j outside the 
context of gyrokinetic theory. The purpose of the present paper is to investigate the higher-order polarization and 
magnetization effects introduced by the guiding-center gyroradius p gc = p 0gc + e# p lgc + • • • within the standard 
format of nonlinear gyrokinetic theory. 

The remainder of the present paper is organized as follows. In Sec. [ill we introduce variational definitions of the 
reduced polarization and magnetization used in modern gyrokinetic theory based on functional derivatives of the 
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gyrocenter Hamiltonian. In Sec. IIII1 we explicitly make use of the higher-order correction p lgc to the guiding-center 
gyroradius p gc to derive explicit expressions for the first-order guiding-center polarization and magnetization. These 
results confirm that the two-step perturbation analysis used in modern gyrokinetic theory [l[ yield a consistent set of 
gyrokinetic Vlasov-Maxwell equations that include all first-order (e^ and ea) polarization and magnetization effects. 



II. GYROCENTER POLARIZATION AND MAGNETIZATION 



The polarization effects in gyrokinetic theory are formally defined in terms of the relation (at a fixed position r at 
time t) between the particle charge density g(r, t) and the gyrocenter charge density p gy (r,i): 

g = g gy - V-P gy , (2) 

where the polarization charge density g po \ = — V • P gy serves as a definition of the gyrocenter polarization P gy . The 
gyrocenter charge density (summation over particle species is implied) 

1? j3— 



g gy = e J F d v (3) 

is defined as a gyrocenter- velocity-space integral of the gyroangle-independent gyrocenter Vlasov distribution F (where 
o = 2ir J dvtt dp with Jacobian J to be defined below and the gyrocenter gyroangle integration has been performed 
explicitly). The gyrocenter polarization 

Pgy = Pgc + 66 Pi gy + • • • , (4) 

on the other hand, is defined in terms of the guiding-center polarization 

P gc = e B Pigc + ■ ■ ■ , (5) 

which vanishes in a uniform magnetized plasma, and the first-order gyrocenter polarization esPi gy , which vanishes in 
the absence of field fluctuations. 

The magnetization effects, on the other hand, are formally defined in terms of the relation (at a fixed position r at 
time t) between the particle current density J(r, t) and the gyrocenter current density J gy (r, i): 



<9P 

J = J g y + -^f + cVxM gy , (6) 



where the gyrocenter polarization current J po i = dP gy /dt is defined in terms of the gyrocenter polarization (j4]) and the 
magnetization current J m ag = c V X M gy serves as a definition of the gyrocenter magnetization M gy . The gyrocenter 
current density 

Jgy = Jgc + £5 Jlgy + • • • (7) 

is defined in terms of the guiding-center current density 

Jgc = e / F -g- d 3 v = J 0gc + e B Jigc + • • ■ , (8) 

which includes first-order guiding-center drifts, while the gyrocenter magnetization 

M gy ee M gc + eg Mi gy + • • • (9) 
is defined in terms of the guiding-center magnetization 

M gc ee M 0gc + e B Migc + ■ ■ ■ , (10) 
where the non-vanishing zeroth-order guiding-center magnetization 



M 0gc = - ( / JlFd 3 v) b (1.1) 

is a well-known result foL [Tl|. 
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Note that the definitions ([2]) and © ensure that the gyrokinetic Vlasov-Maxwell equations satisfy the gyrocenter 
charge conservation law dg gy /dt + V • J gy = 0, since polarization effects dg po \/dt + V • J po i = conserve charge iden- 
tically while the magnetization current density is divergenceless V • J mag = 0. We also note that the quasineutrality 
condition g = used in gyrokinetic theory [l[ can now be expressed (up to order es and es ) as 

g gy - V-P gc = e s V-Pigy, (12) 

where the left side contains the standard polarization effects associated with the lowest-order guiding-center gyroradius 
Pogc as wei l as the higher-order polarization effects associated with the first-order guiding-center gyroradius p lgc and 
gradients Vp 0gc . 



A. Functional Definitions 



The gyrocenter polarization (Q} and magnetization © can be defined as functionals of the gyrocenter Vlasov 
distribution F and functional derivatives of the gyrocenter Hamiltonian as follows [l], fl2| - Here, the gyrocenter 
Hamiltonian 

H gy = H ac + e <5 H igy + e S H 2gy + •" ( 13 ) 

is formally expressed as an asymptotic expansion in powers of €b and eg- The unperturbed gyrocenter Hamiltonian 
is defined as the guiding-center Hamiltonian [4| 

H gc = jvl + TlB, (14) 

where higher-order corrections (e^, for n > 1) can be made to vanish Q. The first-order gyrocenter Hamiltonian (in 
the Hamiltonian representation of gyrokinetic theory in which the guiding-center Poisson bracket is left unperturbed) 
is defined as [l3| 

tfi g y - (e^(X + Pgc ,t) - ^(X + p^t) ■ (^ + ^)), (15) 

where the first-order electromagnetic potentials ($i, Ai) represent the electromagnetic-field fluctuations that perturb 
the guiding-center plasma, and the guiding-center (unperturbed) evolution operator d gc /dt = d/dt + { , H gc } gc is 
defined in terms of the guiding-center Hamiltonian (fT4"|) and the guiding-center Poisson bracket { , } gc (to be defined 
below). The exact expression for the second-order gyrocenter Hamiltonian H2 gy has been given elsewhere [l|, [H[ and 
will not be needed in what follows since we will not be concerned with first-order (es) gyrocenter polarization and 
magnetization. 

The first-order gyrocenter Hamiltonian (|15|) , which can be expanded as H\ gy = H^ y + es -^i g y H , contains the 

lowest-order terms [lH 

K0) - - - c A lsc .(vf b + n d -f d 



ee (e$ lgc - ^A lgc - (*j ||b + ft^p)), (16) 



where 4>i gc = $i(X + p Qgc ,t) and Ai gc = Ai(X + p 0gc ,t), as well as the first-order guiding-center corrections 



rr(l) = __ / _ 
a lgy - \Plgo 



e 



dp, 



e E lgc + - I U|| b X B 1±gc + VL x B 1[]sc b 

c \ ou 

(i)y 



e / . / d gc X d gc p 0g , 



c 



A *-i = ^- + - s 5 at ; < i7) 



where Ei gc = — V$i gc — c 1 dAi gc /dt and Bi gc = V X Ai gc = -Bi|| gc b + Bi^ gc , while dg^X/dt and d^J p 0gc /dt 
denote first-order (es) corrections to v\\ b and Fldp 0gc /d8, respectively. 
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Gyrocenter polarization 



The gyrocenter charge density and the gyrocenter polarization are defined as [12J 



V ■ Pev = 67 1 



F 



where guiding-center polarization effects ([5]) arc defined by the identity 



V-P g c = 



F 



(<5 3 (X + p gc -r)) - ed 3 (X-r) 



d b z. 



(18) 



(19) 



By expanding Eq. (|19|) in powers of p gc and integrating by parts to eliminate i5 3 (X — r), we obtain the guiding-center 
polarization 



gc 



F <p_) d A v - V 



F (pgcPgc) d 6 v) + ... 



(20) 



where dipole ((p gc )) and quadrupole ((p gc p gc )) contributions are shown. By using the fact that (p 0gc ) = 0j i- e -> 
Eq. (|20p has a vanishing zeroth-order term, the first-order guiding-center polarization is defined in terms of the 
functionals 



Pi, 



F 



(Pigc) - V - (( 



PdgcPogc 



POgcPogc 



■VF 



}d 3 v, 



(21) 



where we have grouped terms that directly involve the background magnetic-field nonuniformity and the term that 
directly involves the spatial gradient of the gyrocenter Vlasov distribution F. The latter term (whose ordering is 
assumed to be comparable to es) is generated by the Taylor expansion of the integral 



F (<5 3 (X + p 0gc -r) ) d 6 z 



F 



JlB 
2mSl 2 



Vi F 



d 3 v, 



(22) 



where p Qtrc is assumed to be spatially uniform. We will return to Eq. (|21[) once we have obtained an expression for 
the first-order guiding-center gyroradius p lac in the next Section. 



2. Gyrocenter magnetization 



The gyrocenter current density and the gyrocenter magnetization are defined as 



<9P S 



dt 



cVxM gy ee e s 



F 



SH, 



gy 



<$Ai(r,t) 



5H 



igy 



5Ai(r,i) 



d 6 z 



where guiding-center magnetization effects (JTOJ) arc defined by the identity 



<9P e 



dt 



c V X M„ 



(ig C X 



<5 3 (X + p g 



d gc p 



gc c3 



<5 d (X + p g 



(23) 



d 6 z, (24) 



with the guiding-center current density J gc defined by Eq. (jHJ . After expanding Eq. (|24|) in powers of p gc and carrying 
out several manipulations (see Appendix A for details), we obtain the guiding-center magnetization 



1 / dgcpgc 



+ (Pgc) X 



dt 



d 3 v, 



(25) 



which is naturally split into an intrinsic contribution associated with (1/2) (p gc X d gc p gc /dt) and a moving-electric- 
dipole contribution (p gc ) X d gc x./dt. Note that the moving-magnetic-dipole contribution to the guiding-center polar- 
ization (|2"U1) is a relativistic effect [lij which falls outside the scope of the present work. 
The guiding-center magnetization has a well-known zeroth-order intrinsic contribution [ll[ 



M 0gc EE — 



2c 



(26) 
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from which we recover the classical result where we used the lowest-order identities dp „ c /dd = p „ c X b and 



lPo s 



2p£?/(mil 2 ). The first-order contribution 



M 



lire 



F 




d 3 v 



F 



(Pigc) x [v\\b 



(27) 



on the other hand, is decomposed in terms of the first-order corrections (the first two terms) to the zeroth-order 
guiding-center intrinsic magnetization (|26l) . while the third term yields the moving-electric-dipole contribution since 
it involves the gyroangle-averaged first-order guiding-center gyroradius (pw) appearing in Eq. (|2Tj) . 



III. GUIDING-CENTER POLARIZATION AND MAGNETIZATION 

In the present Section, we make use of the higher-order corrections to the guiding-center transformation 0, 0] for the 
purpose of determining the first-order guiding-center gyroradius p lgc used in the first-order guiding-center polarization 
(f2~Tj) and the first-order guiding-center magnetization (|27l) . We use Northrop's macroscopic interpretation of the 
small parameter which, for finite macroscopic length Lb, allows us to use e _1 ~ e as an ordering parameter. Hence, 
the guiding-center dynamical reduction is generated by the near-identity phase-space transformation 

(1 f)C a \ 
G a 2 + £ -g^J + (28) 

where e ~ e _1 is used an ordering parameter and the phase-space vector field G„ is said to generate the phase-space 
transformation at order n > 1. 



A. Guiding-center phase-space transformation 



The guiding-center coordinates z a = (X, 0) are defined (up to first order in es) as Q 



X = 

v u = 
A* = 



x - ep 

v 



o Po v Po o 7T~ + 



/'n 



1 - ep Q - [b-\7b 

1 - e p || (t + a 1 : Vb 

e p • R + e po,! fa 2 : Vb) 



p 5 



<9p 

t + ai : Vb) , 



6 Po • I Mo V In S - 



de 



o|| 



2p .B 



b-Vb 



b-Vb 



(29) 
(30) 
(31) 

(32) 



where p || = u o||/( TO ^) denotes the "parallel" gyroradius, r = b- V X b denotes the background magnetic torsion, R 
denotes the gyrogauge vector field, and expressions for the gyroangle-dependent dyadic tensors (ai, a 2 ) are not needed 
in what follows. Here, the second-order spatial component GJ is expressed as 



G 2 — G 2 || b 



P\\ T Po + 5 ( G i - MPo-VlnB 
where 

G 2 |, = b-Gf = 2 p || (b-Vb 
The first-order components (G^Gf), defined in Eqs. ([3T 



9p 1 / r e , _ t>A 9 Po 
— + -[G, +Po .Rj_. 



po B 



a 2 : Vb 



and used in Eqs. (J2HJ) and (j3"3")l . are 



Gt = po- po VlnB 



B 



b-Vb 



G? 



Po • R + Po|| (a 2 : Vb) 



dPo 

d9 



Mo Poll 



VlnB 



ai : Vb , 



2p B 



b-Vb 



(33) 
(34) 

(35) 
(36) 
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The remaining first-order component defined in Eq. (|30[) 

G? = - wbn Po ■ (b- Vb) + ^ (r + ai : Vb) (37) 

ensures that the first-order correction 

H lgc = -G?-q£t = MPo-VB - mv\\ G? - B _? = (38) 

to the guiding-center Hamiltonian (|14[) vanishes identically. The same construction algorithm can be applied at higher 
order (i.e., H ngc = — G"dH sc /dz a = for n > 1), so that the simplicity of the guiding-center Hamiltonian can be 
preserved in the form of Eq. (fT4f to all orders in e Q . 

The Jacobian for the guiding-center phase-space transformation (|2"9")) - (f3"2"l) is constructed from the Jacobian for the 
local particle phase-space coordinates J7q = B /m according to the formula 



d B 

J = Jo — e tt— {Jo -?) + • ■ • = — 

oz a m 



. B \ B fdG? dGt dG{ 

V • — Po -~— + -fr 1 - + 

~~\ J m \ ovu on do 



R B* 

= _(i + ep T + ...) = JL. (39) 
m m 

Hence, the near-identity guiding-center phase-space transformation is noncanonical since Bjj ^= B. 

As a result of the guiding-center transformation (I2"!)l) - (f3"2")) . the guiding-center phase-space Lagrangian is expressed 

as 

/ g B \ B 

r gc = — A + mv\\ b - e /i — R* • dX + e /u — d0 - # gc df, (40) 

\ 6 C i ' / i ' 

where R* = R + (r/2) b and the guiding-center Hamiltonian H gc is given by Eq. (| 14[) . The guiding-center Poisson 
bracket { , } gc constructed from the symplectic part of the guiding-center phase-space Lagrangian (39 is i 

r- ittfdFdG dF dG\ B* /_* _<9G dF _„\ cb _„„ 

{F, G} gc - e -i _ _______ _ + __ . V*F — V*G - e -— • V*F x V*G, 41 

B \80 d[i d^i 88 J mB?, \ dv\\ dv\\ J eB^ 



where V* = V + R* d/89, 



B* = B + B ( e p|| V X b - e 2 £^ V X R* + • • • ) , (42) 



and _j| = b-B*. Note that, under the gyrogauge transformation 6 ^ 9' = 9 + ^>(X), the vector R is gyrogauge- 

dcpcndcnt (i.e., R > R' = R + Vt/>), while the curl of R is gyrogauge-invariant (i.e., V X R' = V X R). Hence, the 
guiding-center phase-space Lagrangian (|40[) and the guiding-center Poisson bracket (|4ip are both gyrogauge-invariant 
[J| since the combinations d9 — R • dX and V + Rd/89 are gyrogauge-invariant. 

Lastly, we can now write expressions for the velocities d gc ~K/dt and d gc p Q /dt to be used in evaluating the first-order 
guiding-center magnetization (|25p . First, the guiding-center velocity 

^ _ {X, H gc } gc = v l{ b + x (/x VS + mujf (b- Vb)) _ «„ b + e B ^ (43) 

is expressed in terms of the zeroth-order motion along a magnetic-field line and first-order guiding-center drifts. 
Second, the gyration velocity 

— - = {p , i/ gc } gc = n w + «i|b-v po + ••• = n + + ••• (44) 

is expressed in terms of the zeroth-order perpendicular particle velocity and its first-order correction. The guiding- 
center velocities ([4"3"f and (|4"4")l appear in the first-order gyrocenter Hamiltonian (|T5|) and, consequently, they contribute 
to the guiding-center magnetization (|25[) . 
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B. First-order corrections to the guiding-center gyroradius 



So far we have not made a distinction between the gyroradius in particle phase space (labeled p) and the gyroradius 
in guiding-center phase space (labeled p KC ), which are respectively defined as 



P = x - T gc X = p + € p l + • 

Pgc = "'"gc X X = Pogc + e Plgc 



(45) 



In Eq. (|45l) , the pull-back operator T gc transforms a function on guiding-center phase space into a function on particle 
phase space, while the push-forward operator T^ 1 transforms a function on particle phase space into a function 
on guiding-center phase space. From these definitions, we therefore obtain the relation between the particle and 
guiding-center gyroradii 



T g~cV 



eG? 



dp 

dz a 



Poz 



(46) 



so that the lowest-order guiding-center gyroradius is p 0gc = p is identical in both phase spaces. The first-order 
gyroradii p 1 and p lgc = p 1 — G" dp Q /dz a are not identical, however, and it is important to use the proper first-order 
gyroradius in order to obtain the correct first-order polarization and magnetization. 

In particle phase space, the guiding-center position X is expressed in terms of the guiding-center transformation 
([25)1 as X = x — e p — e 2 p 1 + • • • = T gc x, where the first-order displacement G* = — p Q defines the lowest-order 
gyroradius, while the second-order correction is 



Pi 



— G? - - 



= (^2^o- VmS - P\\ T 
The gyroangle-averaged particle gyroradius (1471) is 



Po 



1 



-p -Vb-p - G 



2|| 



(Pl) = 



fiB 
2mVL 2 



V-b b + V ± lnB 



= - V- 



fiB(b- Vb) 
2 mil 2 



(47) 



(48) 



where we used (G 2 \\) = from Eq. (l34l . 

In guiding-center phase space, on the other hand, where the particle position x is expressed in terms of the inverse 
guiding-center transformation as x = X + p 0gc 



Pi 



gc 



= —Gn 



Po' v Po - 



Pi g , 
Gl 



Tg^X, the first-order correction is 



fi dp 



G'[ 



dPa 
d6 



n Tt\ dpQ 



where the same generating vector fields (Gi, G2, ■ • • 
Po' v Po = - (/xpo'VlnB 



'II 



dpp 



G 



'2\\ 



- p -Vb-p , 



dji 

G^ + 2/f/' 

• ) are used and we inserted the identity 

p -Vb-p ) b 



(49) 



^-(po-R)%- 



89 



in obtaining the last expression for p lgc . 

Lastly, we note that the first-order particle and guiding-center gyroradii ([4"T]) and (14T)1) are both gyrogauge-invariant. 
Moreover, according to Eqs. (|21l) and (|27|) . it is the guiding-center first-order gyroradius vector (|49|) that must be used 
and, hence, we will use (Pi gc ) to compute the first-order guiding-center polarization ([21]) and p lgc = p lgc — (Pi gc ) to 
compute the first-order guiding-center magnetization ([27[) . 



C. First-order Guiding-center Polarization 

The gyroangle-averaged guiding-center gyroradius (j49[) is 
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where we used Eq. (|48|) . When we combine these results into the first-order guiding-center polarization (|2T|) , we obtain 



1 



cVx [- 



-Jib Fcfu 



(51) 



This expression combines the classical first-order (dipole) guiding-center polarization obtained previously Q as well 
as the first-order quadrupole guiding-center polarization. 



D. First-order Guiding-center Magnetization 



The calculation of the first-order guiding-center magnetization (|27l) requires the guiding-center push-forward of the 
particle velocity v = dx/dt, which is defined as 



V gc 



T gc rv 



cL c X i d gc p 



where 



V 



dt 



Ogc 



dt 



"gc FOgc 



= V 



Ogc 



e V 



lgc 



vub + n 



dp 



Ogc 



dt dt " ' 89 

denotes the lowest-order guiding-center (particle) velocity and its first-order correction is 

Vigc = (V lgc ) + v lgc , 

where the gyroangle-independent part is 



and the gyroangle-dependent part is 



"gc Pi g 

dt~ 



dt 



"gc KOgc 

dt 



dt 



x b 



n 



dp 



lgc 



06 



(52) 



(53) 



(54) 



(55) 



Here, we easily show that the guiding-center velocity V gc satisfies the following identities: (V gc ) = d sc X./dt and 
H gc = (m/2) |V gc | 2 . 

Lastly, by inserting Eqs. (|49p -(|50 |) and (|55|) into Eq. fj 2 7|> . we obtain the first-order guiding-center magnetization 
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Mi. 



P\\ F 



H r b 



BVx B~ b 



B 



bx b-Vb 



d 3 v. 



(56) 



We note, however, that the first-order guiding-center magnetization (|56[) vanishes if the gyrocenter Vlasov distribution 
F is a Maxwellian in vu . In most applications of gyrokinetic theory, the zeroth-order guiding-center magnetization 
(|2"6]) is therefore sufficient. 



IV. SUMMARY 



We have presented the derivation of guiding-center polarization and magnetization effects as a simple extension of 
the standard form of modern gyrokinetic theory by keeping higher-order corrections to the guiding-center gyroradius 
(14T)1) in the first-order gyrocenter Hamiltonian (|T5|) [e.g., Eq. (TTT)) ]. 

From the variational derivations of the guiding-center polarization and magnetization, we have recovered the classical 
first-order guiding-center polarization (|51[) and the classical zeroth-order guiding-center magnetization (|26p . These 
results confirm that the two-step perturbation analysis used in modern gyrokinetic theory yield a consistent set of 
gyrokinetic Vlasov-Maxwell equations that include all first-order polarization and magnetization effects. 
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Appendix A: Reduced Polarization and Magnetization 



In this Appendix, we present a brief summary of the derivation of reduced polarization and magnetization effects 
fill induced by a general near-identity phase-space transformation z — > z = %z. The dynamical reduction introduced 
by this phase-space transformation yields the reduced Vlasov equation 



dt j dz a \ dt 1 ' 



(Al) 



where J denotes the Jacobian of the transformation Te and the reduced Hamilton equations d e ~z a /dt = {z a , H} e are 
represented in terms of a reduced Hamiltonian H and a reduced Poisson bracket { , } e (not necessarily canonical). 

Reduced polarization and magnetization effects are associated with the reduced displacement p e = T ( T 1 x — x. We 
begin with the push- forward derivation of the reduced polarization generated by the phase-space transformation T e . 
The particle charge density (summation over particle species is implied) 

(A2) 



(A3) 



and the reduced polarization charge density £ pol = — V - P, where the reduced polarization is defined as a multipole 
expansion (dipole + quadrupole + ■•■) associated with increasing powers of ~p t : 



g = e / / 5 d (x + p e - r) d b z = g - V-P 
is expressed in terms of the reduced charge density 

q = e [ J8 3 (x~r)d 6 z = e [ J d 3 v, 



f Pe d 3 V 



V- I - / fp e p e d"v 



(A4) 



Next, we consider the push-forward derivation of the reduced magnetization generated by the phase-space trans- 
formation Te. This derivation uses the push-forward transformation of the particle velocity v = dx/dt: 



T7 x v 



Tr x x 



d e x 
~dT 



The particle current density 



dePe 



dt dt 

is thus expressed in terms of the reduced current density 



<5 3 (x + p e -r) d 6 z = J 



dePe 

dt 



dP 

~dt 



V X M, 



d e X 3 _ 



(A5) 



(A6) 



(A7) 



where J 



pol 



dP/dt is the reduced polarization current density, and the reduced magnetization current density 



i V X M. If we expand Eq. in powers of p e , we obtain 



/ 



dt 



dePe 

dt 



d 3 V 



fPe 



dt 



dePe 

dt 



d 3 V 



(A8) 



Next, we take the partial time derivative of the reduced polarization (|A4[) . with the reduced Vlasov (|Aip . we can 
insert 
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into Eq. (|A8I) and obtain 



J = J+ m 



f 



dT p ^ ~ p * -dt 



1 ( deP e _ _ d e p € 



dt 



dt 



(A9) 



Lastly, by using the identity V • (BA — AB) = V X (A X B), for two arbitrary vector fields A and B, we obtain an 
expression for the reduced magnetization 



M 



/ 



Pe X 



1 d e p e 

2 dt 



~~dt 



d 3 v 



(A10) 



In Eq. (|A10p , the term associated with (1/2) p € X d e p e /dt represents the intrinsic reduced magnetization while the 
term ~p e X d e x./dt represents the moving-electric-dipole contribution. 
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